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Abstract 

We consider the insurance company as a physical system which is immersed in its environment 
(the financial market). The insurer company interacts with the market by exchanging the money 
through the payments for loss claims and receiving the premium. Here in the equilibrium state 
we obtain the premium by using the canonical ensemble theory, and compare it with the Esscher 
principle, the well known formula in actuary for premium calculation. We simulate the case of car 
insurance for quantitative comparison. 
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I. INTRODUCTION 



From the first of nineteen century the economists have tried to find a way to use the 
formalism of physical theories in mathematical modeling of economy. The concepts of utility 
function and economical equilibrium entered into the financial theories in correspondence 
with potential energy and mechanical stability [lj. The physicists have paid less attention to 
this subject til. last decade Q. Recett.y they have more regards for studying the dyuauucs 
of the stock markets |2J, |3J. The main source of their interest comes from fluctuation of stock 
prices with the time. Understanding this behavior enable us to manage the strategies for 
trading 4]. The statistical mechanics has appeared as a powerful method for exploring the 
price dynamics [l|, E , 0, Q ■ 

The insurance is one of the important parts of the financial market with respect to trading 
risks. In the specified period of time, the results of a loss events (risks) that occur for the 
insurant are covered by insurance company and insurant also pays an amount of money 



to the insurer, which is called premium p, |f|. In the competitive financial market the 
calculation of premium is very complicated. The premium is affected by random nature of 
risks and also variation in number of the insurants. The latter may be decreased due to an 
increment in premium and is increased by reduction of premium. The actuary is a branch of 
mathematics that studies the relation between paid and received money in order to assign 
premium to a category of risks. 

Nowadays the physicists also pay attention to this part of the financial market. They 
look for a new way for the premium calculation on the basis of physical concepts j?l 0| • 

In this paper we are intended to explain how the equilibrium statistical mechanics may 
be used for premium forecasting. In the next section we describe the analogy between the 
financial market and physical system in contact to heat reservoir. A method for equilibrium 
insurance pricing will be suggested based on the canonical ensemble theory. In the last 
section the special case of car insurance is simulated and we compare our method with 
Esscher principle which is the famous method in actuary and asset pricing jlJEli, U, 
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II. DESCRIPTION OF THE MODEL 



The financial market is combination of large number of economic agents which are inter- 
acting with each other through buying and selling. We consider the behavior of one of the 
agents for example an insurance company; all other agents may be regarded as its environ- 
ment. The agent exchanges money when interacts with its environment. We suppose the 
financial market is a closed system, the environment will absorb the money that the agent 
loses and supplies the agent's incomes. In equilibrium state a parameter that we denote it 
by (3 has the same value for the agent and its environment jisl ]. 



P Agent P Environment ( 1 ) 

This parameter should be related to macroscopic properties of an economic agent. Later 
on, we signify it for an insurance company in terms of its initial wealth, the mean claim size 
and the ruin probability. 

Like a physical system in contact with a heat reservoir we can nominate the following 
probability for E(t), the net money w 
at the end of specified time interval 



rich is exchanged between an agent and its environment 

3, 



-f3E(t) 

?rm)) = ^z^y (2) 

The sum is over all possible amounts of exchanged money in equilibrium state. We assume 
here that the agent's money is considerably lesser than the money in the environment. 
The above result is confirmed by simulation jl^l Q] and is also inferred from empirical 



data 
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For the special case of car insurance in the time interval t, the insurance company receives 
premium from I(t) insurants and covers the loss results for N(t) car accidents. In each 
accident insurant charges the company for Xj amount of money. Thus the exchanged money 
at the end of interval is, 

N(t) 

E{t) = pl(t) - £ X, (3) 

Where p is the premium. In the insurance terminology we call this money the surplus of 
company. The /(£), N(t) and Xj are random variables which their probability distribution 
should be inferred from empirical data of the insurance company. 



Premium calculation principle is a rule that assign to any distribution function (corre- 



spond with loss events) a real number 



5] 



Biihlmann in his brilliant papers on economic 



premium principle stated that premium calculation principle should also depend on market 



conditions in addition to loss distribution 



la. Il9| . He constructed a model for whole insur- 



ance market with many agents; like insurance companies, reinsurance companies and buyers 
of insurance. He concluded that the sum of all incurred losses for all agents in the market 
shows the effect of surrounding market and it should be included in premium calculation 
principle. 

In this work we introduce a new way for premium calculation. The number of insurants 
enters in our model as a reminiscence of the market conditions, the competition in the 
market changes it randomly as we mentioned above. 

The insurance company avoids of financial failure. Premium should be calculated so that 
the surplus of company at the end of specified interval should be zero at least. This condition 
may be satisfied in mean. 

<E(t)>=J2Pr(E(t))E(t) = 0. (4) 

In real case premium is greater than the value obtained from eq.|3]because administration 
expenditures must be compensated and company also likes to profit, anyhow we aren't 
considering these matters. 

If the number of insurants becomes constant then by simple calculation eq. |U reduces to 
what is known as the Esscher principle. 

p J o = • ( 5 ) 

Where Z(t) is the aggregate loss in the time interval t, it is defined as, 

N(t) 

Z(t) = £ X j (6) 

3=1 

Unlike to the net premium that is obtained by averaging the aggregate loss, in the Esscher 
principle the premium is calculated by exponential tilting of loss events. The parameter (3 
appears to scale the weight of each event. There is no unique way for computing the (3 
parameter in Esscher principle; some authors use the optimization procedure |5| and others 
relate it to the risk tolerance of the insurance company [3]. The way that we adopt in this 
paper for determination of this parameter is also applicable here. 



The Esscher principle as an efficient method for premium computing appears as a special 
case of our method. The difference between it and our method is shown in the next section 
numerically. 



III. SIMULATION RESULTS 

To simulate our model, it is necessary to specify the distributions for the random quan- 
tities which are appeared in expression for the surplus, eq. El 

By assumption the number of loss claims by insurants, N(t) has Poisson distribution with 
the mean Xt. It depends only on the length of interval t and also is statistically independent 
of the number of claims in any previous interval. 

Pr{N{t + s)-N{s)=n) = ^y (7) 

71/. 

The important property of Poisson distribution is that the times between successive 
claims are independent and identically exponentially distributed (the Erlang's model) with 
the mean 1/A. It is what we encounter in most of the real cases. 

The model for the claim size is constructed based on the knowledge and experience of the 
insurance company in addition to data from the past. The model should provide a balance 
between simplicity and conformity to the available data. 

The purpose of this paper is explanation of our method; in this respect for simplicity we 
restrict ourselves to the exponential distribution for the claim size. Otherwise we are forced 
to do a lengthy numerical computation that don't correspond to our goal. 

Pr(X = = -e-*/" (8) 

The n is the mean of the claim size. The product of the Xt and fi is the net premium, 
Po. The insurance company which offer this value or less for the premium will be busted 
ultimately 

The distribution for the number of insurants has not been modeled since now because it 
doesn't be involved in the traditional methods for calculating the premium. The eq. |U shows 
its importance in our method. However we assume the number of insurants is distributed 
uniformly with the mean, J , and standard deviation, %10I . This assumption based on 
empirical data of the Iran Insurance Company 



so 
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FIG. 1: The annual profit or surplus after 300 days versus the premium for large (3. The premium 
and surplus are measured in terms of net premium. 
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FIG. 2: The premium as a function of fit. The results for t = 10, 100 and 1000 is plotted. 

All the necessary parameters for tuning the mentioned distribution functions are obtained 
: rom statistical reports of Iran Insurance Company for the case of car insurance in year 2000 
7 or independence of the results to any monetary units, the initial wealth and premium 
are expressed in terms of the mean claim size and net premium. The following results are 
averaged over 10000 different possible states of the modeled insurance company. 

The variation of annual profit or surplus after 300 days with premium for large [3, is 
shown in the fig. Q The desired value for the premium is obtained from calculation of the 
intersection point of approximated curve with the premium axis. For other values of the 
parameter (3 we can apply the same method for computing the premium. 

The premium is also function of parameter (3. In the case of [3 — it is equal to the 
net premium. An increment in (3 increases the premium but for large (3 it approaches to a 
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FIG. 3: The premium versus duration of the insurance contract for the large f3. The squares are 
obtained by simulation and the curve shows the fitting result. 
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FIG. 4: The parameter ln(/3i) as a function of filnell 1 which is introduced in the text. The results 
for U = 1, 10 and 1000 times of the mean claim size are plotted. 

limiting value. The fig. El shows this behavior for three value of the time interval. 

The fig. 01 demonstrates dependency of the premium on duration of insurance contract 
for the large (3. This behavior is expected in the insurance market p. The premium for 
other values of the parameter (3 shows the same behavior in the time. 

The parameter /3 has been ambiguous since now. At this moment we proceed to clear 
its meaning. It can be understood by dimensional analysis of the eqs. 01 and |U that the 
parameter (3 is proportional to the inverse of time. The fig. El also confirms this statement, 
the behavior of the premium as a function of j3t is the same for different time intervals in 
exception of their limiting values. 

The most important quantities which are used by the insurers for their financial decisions 
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are the initial wealth U and the ultimate ruin probability e. The combination of them 
with the mean claim size /i is used in all methods for calculating the premium ^, |(|. The 
parameter (3 is also function of these quantities. The fig H] shows the In (3t as a function of 
the dimensionless parameter —\x In ell" 1 and the initial wealth. This consequence is issued 
from relation of the ultimate ruin probability on premium p, 

£ = — exp 1 . 9 

p up 

And the result for dependency of premium on (3t which is also shown in fig. |21 

The ultimate ruin probability for exponential distribution of claim size can be calculated 
analytically as is seen in the eq. El but in other cases it would be computed numerically 
from the following relation [6j], 

e = l- V(l -^)(^) fc F e * fe (f/). (10) 
fc =o P P 

Where F e (x) is related to the claim size probability F(x), 

F e (x) = - [ X [l-F(y)]dy (11) 
H Jo 

And F* k (U) is the /c-fold convolution of F e (U) with itself. 

The fig. |U also indicates if the initial wealth becomes large for a fixed time interval the 
parameter (3 approaches to zero. This means that a wealthy company buys the risks with the 
less premium in comparison to the small companies which offer more expensive insurance. 
This fact conforms to our experience in the financial market. 

As mentioned before the Esscher principle is concluded from our method. The variation 
in number of insurants changes the premium certainly. The fig. displays the difference 
between two methods numerically. 



IV. SUMMARY 



We construct a simple model for the cash flow in an insurance company for the category 
of non-life insurance. A new method is suggested for calculation of premium on the basis 
of the canonical ensemble theory. This method considers the effect of the environment (or 
market conditions) by introducing a new random variable, the number of insurants, into 
the equation of cash flow. In this respect the empirical data of an insurance company is 

8 



0.07 
0.06 
0.05 
V 0.04 

o 

fi, 0.03 
0.02 
0.01 
0.00 

FIG. 5: The premium versus the parameter [3t for 100 days interval. The results of our method 
(canonical ensemble theory) and Esscher principle are plotted for comparison. 

sufficient for take into consideration of environment conditions. There is no need to acquire 
any information about other agents in the market unlike the Buhlmann method. Another 
advantage of this method is ability for extension in the small market case. The latter is 
under investigation by the author. 

The financial backbone of insurer company is concealed in the parameter f3 which play 
essential role in our method for computing the premium. When we omit the effect of market 
conditions by considering constant value for number of insurants the Esscher principle is 
emerged as a result. Difference between our method and the Esscher principle that is 
apparently due to the effects of environment is shown numerically. 

This work has been supported by the Zanjan university research program on Non-Life 
Insurance Pricing No 8243. 
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